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Abstract 

In this paper we study the Geodesic Deviation Equation (GDE) in metric f(R) gravity. We start 
giving a brief introduction of the GDE in General Relativity in the case of the standard cosmology. 
Next we generalize the GDE for metric f(R) gravity using again the FLRW metric. A generalization 
r—i of the Mattig relation is also obtained. Finally we give and equivalent expression to the Dyer-Roeder 

equation in General Relativity in the context of f(R) gravity. 
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In the context of General Relativity (GR), the curvature and geometry of the space-time play a fun- 
(3JT) damental role: instead of forces acting in the Newtonian theory we have the space-time curved by the 

matter fields [1]. The space-time could be represented by a pair (Ai,g) with M. a 4-dimcnsional mani- 
fold and g a metric tensor on it, the curvature is described by the Riemann tensor R, and the Einstein 
field equations tell us how this curvature depends on the matter sources. We can see the effects of the 
^\ curvature in a space-time through the Geodesic Deviation Equation (GDE) [2]- [4]; this equation give us 

the relative acceleration of two neighbor geodesies [5] , [G] . The GDE gives an elegant description of the 
structure of a space-time, and all the important relations (Raychaudhuri equation, Mattig Relation, etc) 
could be obtained solving the GDE for timelike, null and spacelike geodesic congruences. 

o 

Even GR is the most widely accepted gravity theory and it has been tested in several field strength 
regimes, is not the only relativistic theory of gravity [7]. In the last decades several generalizations of the 
• • Einstein field equations have been proposed [8], [9]. Within these extended theories of gravity nowadays a 

subclass, known as f(R) theories, are an alternative for classical problems, as the accelerated expansion 
of the universe, instead of Dark Energy and Quintessence models [10]- [16]. f(R) theories of gravity are 
basically extensions of the usual Einstein-Hilbert action in GR with an arbitrary function of the Ricci 
scalar R [17]-[19]. 



1 Introduction 



In the context of Palatini f(R) gravity the GDE had been studied in [20], where they get the modi- 
fied Raychaudhuri equation. Our aim in this paper is to obtain the GDE in the metric context of f(R) 
gravity, and study some particular cases, we also obtain the generalized Mattig relation which is an useful 
equation to measure cosmological distances. Through this paper we use the sign convention (— , +, +, +) 
and geometrical units with c = 1. 
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2 Field equations in f(R) gravity 



As we mentioned above the modified theories of gravity have been studied in order to explain the accel- 
erated expansion of the universe, among other problems in gravity. A family of these theories is modified 
f(R) gravity, which consists in a generalization of the Einstein-Hilbert action (Lagrangian C[g] = R—2A), 
with an arbitrary function of the Ricci scalar (£[g] = f(R)) [11], [19]. Again, we consider the space-time 
as a pair (M , g) with g a Lorentzian metric on M. ; the relation of the Ricci scalar and the metric tensor 
is given assuming a Levi-Civita connection on the manifold, i.e. a Christoffel symbol. Then the action 
could be written with a boundary term as [21] 

(2.1) 



/ d 4 x^f(R) + 2 f d 3 yey/\h\f'(R)K + S M , 

JV JdV 



here h is the determinant of the induced metric, K is the trace of the extrinsic curvature on the boundary 
dV, e is equal to +1 if dV is timelike and —1 if dV is spacelike (it is assumed that dV is nowhere null), 
and K — 8nG. Finally Sm is the action for the matter fields. Variation of this action with respect to g a ^ 
gives the following field equations [21], which are equivalent for those without boundaries given in [22] 



/ (R) R a ^ 



f(R) 



g a0 + g aP Uf{R) - V a Vpf'(R) = nT a p, 



where □ = V CT V CT , f'{R) = df(R)/dR; the energy-momentum tensor is defined by 



T, 



A I 



a0 



dg 



■ 1V1 ~ 

+ l~, M gaP 



2 SS 



9 Sg 



a/3 ' 



(2.2) 



(2.3) 



being Cm the lagrangian for all the matter fields, we also have the conservation equation V Q ,T Q ' 3 = 0. 
Contracting with g a " we have for the trace of the field equations 



f'(R)R-2f(R) + 3af(R) = KT, 
with T = g a ^T a p. From equation (2.2) we can write 



Rat 



f'(R) 



KT aP + g ali - g a0 nf(R) + V a V fj f{R) 



and from (2.4) 



f'{R) 



R = — — kT + 2f(R) - 3Df'(R) 



(2.4) 



(2.5) 



(2.6) 



It is possible to write the field equations in f{R) gravity, in the form of Einstein equations with an 
effective energy-momentum tensor [15] 



G a R = R a t 



1 



:Rg a p 



KT aP lf(R) ~ Rf'(R)] M [V Q V^/'(i?) - g aP af'(R)} 

9a/3 „ „./ ^ r 



f'(R) 



2f'(R) 



f'{R) 



or 



with 



Gat 



T 



eff - 



f(R) 

[f(R) Rf'(R)] 



(T, 



a/3 



rpeff\ 



9ufs + [V a V^ - g aP U]f{R) 



(2.7) 
(2.8) 

(2.9) 



which could be interpreted as an fluid composed by curvature terms. 
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3 Geodesic Deviation Equation 



Now we give a brief discussion about the Geodesic Deviation Equation (GDE) following [4]- [6]. Let be 
7o and 71 two neighbor geodesies with an affine parameter v. We introduce between the two geodesies 
a entire family of interpolating geodesies s, and collectively describe these geodesies with x a (v : s), figure 
1.1. The vector field V a = is tangent to the geodesic. The family s has rj a = like it's tangent 
vector field. Thus, the acceleration for this vector field is given by [5], [6] 




Figure 1.1 Geodesic Deviation. 

Thus, the acceleration for this vector field is given by [5], [6] 



2„Ot 



~Eh7 2 



(3.1) 



which is the Geodesic Deviation Equation (GDE). Here correspond to the covariant derivative a long 
the curve. We want to relate the geometrical properties of the space-time (Riemann and Ricci tensors) 
with the matter fields through field equations. For this we write the Riemann tensor as [5], [23] 

1 R 

RafijS — CafSjS + ^ (da^Rsp ™ 9a8R~jf3 + gpsR'ja ~ 90fRsa) ~ (9aj9Sp ~ 9a89~tp) , (3-2) 

with Ca/3-yS the Weyl tensor. In the case of standard cosmology (Friedman-Lamaitre-Robertson- Walker 
universe, hereafter FLRW universe) we have the line element 



ds 2 



-dt 2 + a 2 (t) 



dr 2 



1 — kr 2 



+ r 2 d9 2 + r 2 sin 2 9dip 2 



(3.3) 



where ait) is the scale factor and k the spatial curvature of the universe. In this case the Weyl tensor 
C a p~jS vanishes, and for the energy momentum tensor we have 



T a p = (p + p)u a Ufs +pg a p, 
being p the energy density and p the pressure, the trace is 

T = 3p - p. 



(3.4) 
(3.5) 
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The standard form of the Einstein field equations in GR (with cosmological constant) is 



then we can write the Ricci scalar R and the Ricci tensor R Q p using (3.4) 

R = k{ p - 3p) + 4A, 

from these expressions the right side of equation (3.1) is written as [4] 



(3.6) 

(3.7) 
(3.8) 

(3.9) 



with e = V a V a and E = —V a u a . This equation is known as Pirani equation [3]. The GDE and some 
solutions for spacelike, timelike and null congruences has been studied in detail in [4], which gives some 
important result concerning cosmological distances also showed in [23]. Our purpose here is to extend 
these results from the modified field equations in metric f(R) gravity. 



4 Geodesic Deviation Equation in f(R) Gravity 

The starting point are the expressions for the Ricci tensor and the Ricci scalar from the field equations 
in f(R) gravity, equations (2.5) and (2.6) respectively. Using the expression (3.2) we can write 



R, 



i 



2f'(R) 

g a ^S0 — g a s'D 1 p+ gp S V ia — g^Uga) f (R) 



K(Tspg aj — T l/3 g aS + T ia gp S — Tsag^) + f{R){g a - i gsp — gasg^p) 
1 



6f(R) 



kT + 2f(R) + 3D f(R) (g aj gs/3- g a sg~/p), 



where we defined the operator 



(4.1) 



(4.2) 



Thus rising the first index in the Riemann tensor and contracting with V^rpV 6 , the right side of the 
GDE could be written as 



R^v-'V 5 = CfoV^V* + 



K (T S0 5° - T lfj 6f + T<*g 0S - T s a g M ) + f(R)(6«g Sfl - Sfg yP ) 



+ ffiVsp - SfV yP + g ps V° - g M V s a )f{R) 



- Ur + 2f(R) + 3Df'(R)) (S«g sp - Sfg^V^V 5 , (4.3) 

In the following section we explicitly show the steps in order to find the GDE in f(R) gravity using 
FLRW metric, our purpose is compare with the results from GR in the limit case f(R) = R — 2A. 
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4.1 Geodesic Deviation Equation for the FLRW universe 

Using the FLRW metric as background we have 

1 



Ra/3 - 77 



f'{R) 



R 



( f(R)\ 
n(p + p)u a Uf} + ( up -\ — \g af} +V af) f'(R) 

K(3p-p)+2f(R)-2Df{R) 



1 



f'(R) 



(4.4) 
(4.5) 



with these expressions the Riemann tensor could be written as 
1 



R, 



2f'(R) 



k(p + p)(usupg a ~ f - u 7 upg a s + u 7 u a gps - u s u a gpry) 



+ (^p + "y + + a f(R)j {gayggp - gaSgjfl) + (gayVsp - gagDjP + gp&'D 1 a - gpyVs a )f(R) 



(4.6) 



If the vector field V a is normalized, it implies V a V a — e, and 
1 



R afil sV p V 



2f'(R) 
*P , f(R) 



K{p + p)(g ai (uf 3 V 13 )' 2 - 2(u f) V )V( a u 1 ) +eu a u 1 ) 



nf(R) (eg ai - V a V 7 ) + {g ai V sp - g aS V lP + g ps V ia - g Pl V Sa )f{R)V p V & 



up + 



rising the first index in the Riemann tensor and contracting with rf 1 
1 



(4.7) 



2/' OR) 



K(p + p)((u p V l3 ) 2 r ] a - {u p V p )V a ( Ul ^) - (u V^)u a (V 7 rf) + eu a u^) 



(4.8) 



with E — — V a u a , rj a u a = rj a V a = [4] it reduces to 
1 



R^sV^V 6 = 



2f(R) 



1 



2f'(R) 

From the FLRW metric the expression for the Ricci scalar is 



R = 6 



a 




k ' 






- + 






= 6 


H 


a 




a?. 







H + 2H' 



rf. (4.9) 



(4.10) 



where we have used the definition for the Hubble parameter H = -. We see that R is only a function of 
time, and only time derivatives in the operators T> a p will not vanishing. The non-vanishing operators are 

nf(R) = -dlf{R) - md f'(R), 

= -f"(R)R - f"'{R)R 2 ~ 3Hf"(R)R, (4.11) 
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V m = -3Hd f'(R), 
= -3H.f(R)R, 



(4.12) 



V l3 = 2Hg ij d f , (R)+g ij d*f(R), 

= 2H 9lJ f"{R)R + gij (f"(R)R + f"(R)R 2 ), 



(4.13) 



being gij the spatial components of the FLRW metric and R = 8qR. With these results is easy to show 
that the the total contribution of the operators in (4.9) is 

{8°V 6p - 6fV^ + gps V« - g M V s a )f{R)V»V 5 rP = e(5Hf"(R)R + f"(R)R + f"'{R)R 2 )if, (4.14) 
then the expression for R^ s V^rfV s reduces to 



2f'{R) 



n(p + p)E 2 + e [up 



w , f(R) 



2Hf"(R)R 



(4.15) 



which is the generalization of the Pirani equation. In the particular case f(R) = R — 2A the previous 
equation reduces to (3.9). 



The GDE in f(R) gravity from equation (3.1) is 



Dv 2 



2f'(R) 



n(p + p)E + e [up 



W /Cg) 

3 3 



2H,f"{R)R\ 



(4.16) 



as we expect the GDE induces only a change in the magnitude of the deviation vector i] a , which also 
occurs in GR. This result is expected from the form of the metric, which describes an homogeneous 
and isotropic universe. For anisotropic universes, like Bianchi I, the GDE also induces a change in the 
direction of the deviation vector, as shown in [24]. 

4.2 GDE for fundamental observers 

In this case we have V a as the four- velocity of the fluid u a . The afhne parameter v matches with the 
proper time of the fundamental observer v — t. Because we have temporal geodesies then e = — 1 and 
also the vector field are normalized E = 1 , thus from (4.15) 



1 



2f'(R) 



2k P f(R) 



- 2Hf"(R)R 



if the deviation vector is r/ a = £e a , isotropy implies 

De a 



= 0, 



and 



Dt 

D 2 T] a _ d 2 £ 
Dt 2 ~ dfi' 

using this result in the GDE (3.1) with (4.17) gives 



1 



d 2 £ 



dt 2 2f'(R) 



2np f(R) 



2Hf"(R)R 



(4.17) 

(4.18) 
(4.19) 

(4.20) 



G 



In particular with t = a(t) we have 



d _ 1 
a ~ J{R) 



6 3 



(4.21) 



This equation could be obtained as a particular case of the generalized Raychaudhuri equation given 
in [25]. Is possible to obtain the standard form of the modified Friedmann equations [19] from this 
Raychaudhuri equation giving 



1 



3f(R) 



(Rf(R)-f(R)) „ • 

Kp H 3HJ (R)R 



and 



2H + 3H Z 



f'{R) 



K P + 2Hf"(R)R 



(f(R)-Rf(R)) 



f"{R)R + f"'(R)R 2 



(4.22) 



(4.23) 



4.3 GDE for nulll vector fields 



Now we consider the GDE for null vector fields past directed. In this case we have V a = k a , k a k a = 0, 
then equation (4.15) reduces to 



1 



2 „a 



2f'(R) 



K{p + p)E 2 rj 



(4.24) 
1, e a u a = 

e a k a = and choosing an aligned base parallel propagated = fc^V pe a = 0, the GDE (4.16) reduces 



that can be interpreted as the Ricci focusing in f(R) gravity. Writing ry Q = r/e", e a e 
to 



d 2 r] 



1 



-k(p + p)E 2 r\. 



(4.25) 



2f'(R) 

In the case of GR discussed in [4], all families of past-directed null geodesies experience focusing, provided 
K(p + p) > 0, and for a fluid with equation of state p — —p (cosmological constant) there is no influence 
in the focusing [4]. From (4.25) the focusing condition for f(R) gravity is 



f'(R) 



> 0. 



(4.26) 



A similar condition over the function f(R) was established in order to avoid the appearance of ghosts 
[10], [26], in which f'(R) > 0. We want to write the equation (4.25) in function of the redshift parameter 
z. For this the differential operators is 

d dz d 
dv dv dz ' 



(4.27) 



In the case of null geodesies 



d 2 dz d 
dv 2 dv dz \ dv 
dv N 
dz 



(1 + *)=°° 
a 



E 



d^_ 

dz 2 dz dz 2 



dz 
l + z 



da 
a 



(4.28) 



(4.29) 
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with a the scale factor, and oq = 1 the present value of the scale factor. Thus for the past directed case 

a 



dz=(l + z)~~dv = (1 + z)-Edv = E H(1 + zf dv, 
adv a 



dv _ 1 

dz ~ E H{l + z) 2 ' 



dv 

dz 2 



1 



(4.30) 
(4.31) 
(4.32) 

(4.33) 

dz dz dv dt H(l + z) dt ' 

the minus sign comes from the condition of past directed geodesic, when z increases, v decreases. We 
also use ^ = Eq(1 + z). Now, from the definition of the Hubble parameter H 



Then we get 



and 



writing ^f- as 



E H(1 + zf 



dH dv dt dH 



1 , .dH 
H^ 1 + Z ^ +2 



1 dH 



dH dad 

H = — = - — = H, 

dt dt a a 



and using the Raychaudhuri equation (4.21) 

1 



H 



f'(R) 



f ^ + Hf(R)R-^ 
6 3 



H 



then 



d 2 v 



dz 2 E H(1 + zf 
Finally, the operator is 



1 + TO(f-^-^"(«)« 



d 2 ^ 
dv 2 



(EH(l + z)Y 



d 2 i] 

dz 2 + (1 + 2) 



3H 2 f(R) \ 3 6 J \ J J 



dr) 

dz 



(4.34) 

(4.35) 
(4.36) 

(4.37) 



and the GDE (4.25) reduces to 



d 2 ?] 3 
dz 2 + (1 + z) 



3H 2 f'(R) \ 3 6 J 



dt] 



n(p + p) 



dz 2H 2 (1 + zff'(R) 



77 = 0, 



(4.38) 



The energy density p and the pressure p considering the contributions from matter an radiation could be 
written in the following way 

k P = 3Hfo mQ (l + zf + 3Hfo r0 {l + zf, up = Hfo r0 (l + zf, (4.39) 

where we have used p m — and p r = \p r - Thus the GDE could be written as 



fl+V(H,R,z)^- + Q(H,R,z)r, = 0, 
dz* dz 



with 



V(H,R,z) 



4fi m0 (l + zf + 4fi r0 (l + zf + 3f'(R)n k0 (l + zf + AVL DE - 

(i + z)(n m0 (i + zf + Mi + zf + f{R)n k0 (i + zf + n DE ) 



(4.40) 



(4.41) 
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Q(H,R,z) 



3fi m0 (l + z) + 4Q r0 (l + zf 

2(n m0 (i + zf + firo(i + zf + + + M ' 



and given by the modified field equations (4.22) 
1 



ff 2 = 



f(R) 



Hfo m0 (l + zf + H$n r0 (l + zf + 



(Rf'(R)-f(R)) 
6 



- Hf"(R)R 



H 2 = H% 



1 

7W) 



(n m0 (i + zf + n r0 (i + zf + n DE ) + n*(i + zf 



where 



and 







DE 



H 2 



(Rf(R)-f(R)) 
6 



Hf"(R)R 



kO 



H 2 n 2 ' 
n a 



(4.42) 



(4.43) 



(4.44) 



(4.45) 



In order to solve (4.40) it is necessary to write R and H in function of the redshift. First we define the 
operator 

d dz da d d .„> 

= -{l + z)H—, (4.46) 



dt da dt dz 



dz 



then the Ricci is [16] 



R = 6 
= 6 
= 6 



2H Z + H 



r dH 



2H 2 -{l + z)H— + k{l + z) 2 



if we want H = H(z) is necessary to fix the form of H(z) or either a specific form of the f(R) function. 
This point has been studied in [16] and the method to fix the form of H(z) and find the form of the f(R) 
function by observations is given in [27]. 

In the particular case f(R) = R — 2A, implies f'(R) = 1, f"(R) — 0. The expression for fljjE re- 
duces to 



Q.DE = 



1 



{R~R + 2A) 
6 



A 



Ha, 



(4.47) 



then the quantity £Ide generalizes the Dark Energy parameter. The Friedmann modified equation (4.43) 
reduces to the well know expression in GR 



h 2 = h 2 [n m0 (i + zf + a rQ (i + zf + n A + n k (i + zf] , 

the expressions V ', and Q reduces to 

4Sl r0 (l + zf + (7/2)f2 ro0 (l + zf + 3O fe0 (l + zf + 20 A 



V(z) 



(i + z)(firo(i + zf + n m0 (i + zf + n fc0 (i + zf + n A ) 



20*) = 



2^ r0 (l + z) 2 + (3/2)» m0 (l + z) 

n r0 (i + zf + n mQ (i + zf + Q fc0 (i + zf + n A ' 



(4.48) 

(4.49) 
(4.50) 
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and the GDE for null vector fields is 

d 2 r) 4q. (l + z) 4 + (7/2)n m0 (l + z) 3 + 30fcp(l + z) 2 + 20 A d?? 
dz 2 (1 + z)(O r0 (l + z) 4 + m0 (l + z) 3 + feo (l + z) 2 + ft A ) dz 

2q. (l + z) 2 + (3/2)O 7i0 (l + z) 



o r0 (i + z) 4 + q n0 (i + z) 3 + n k0 (i + z) 2 + n A 



rj = 0. (4.51) 



The Mattig relation in GR is obtained in the case Q\ = and writing Slfco = 1 — ^mo — ^ro which gives 
[4] 

fv 6 + n m0 (i + 7z) + n r0 (i + 8z + 4z 2 ) dn 3n m0 + 4n r0 (i + z) =Q 

dz 2 2(1 + z)(l + Q m0 z + tt r0 z(2 + z)) dz 2(1 + z)(l + Vt mQ z + f2 r0 z(2 + z)) ^ ' ' 
then, the equation (4.40) give us a generalization of the Mattig relation in f(R) gravity. 

In a spherically symmetric space-time, like in FLRW universe, the magnitude of the deviation vector 
■q is related with the proper area dA of a source in a redshift z by dr\ oc %/ dA, and from this, the definition 
of the angular diametral distance d A could be written as [28] 

[dA 

with dfl the solid angle. Thus the GDE in terms of the angular diametral distance is 
d 2 d m) ^ 4r> m0 (l + z) 3 +4a, (l + z) 4 + 3/'(i?)^o(l + ^ 2 +4^- dd fW 



dz 2 (i + z)(n m o(i + z) 3 + n rQ (i + z)4 + f>{R)n ko {i + z) 2 + n DE ) dz 

, 3q n o(l + z)+4» r0 (l + z) 2 m) = ^ 

2(ano(i + ^ 3 + a.o(i + 2) 4 + /'(i?)^o(i + ^) 2 + ^D B ) A ' l ' j 

f ( J}\ 

where we denote the de angular diametral distance by d A , to emphasize that any solution of the 
previous equation needs a specific form of the f(R) function, or either a form of H(z). This equation 
satisfies the initial conditions (for z > zq) 



d f A iR) (z,z ) 



= 0, (4.55) 

z=z 



dd 



f(R) 

A 



dz 



-(Z,Z ) 



Ho 



H{zv){\ + z Q y 



(4.56) 



being H{zq) the modified Friedmann equation (4.43) evaluated at z = zq. 

5 Dyer-Roeder like Equation in f(R) Gravity 

Finally we get an important relation that is a tool to study cosmological distances also in inhomogeneous 
universes. The Dyer-Roeder equation gives a differential equation for the diametral angular distance d A 
as a function of the redshift z [29] . The standard form of the Dyer-Roeder equation in GR can be given 
by [30], [31] 

(1 + z) 2 J-(z)^ + (1 + z)g{z)^ + H(z)d A = (5.1) 

with 

T(z)=H 2 (z) (5.2) 
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JIT 

G(z) = (l + z)H(z)^+2H 2 {z) (5.3) 
dz 

KM = ^a„o(l + zf (5.4) 

with a(z) is the smoothness parameter, which gives the character of inhomogeneities in the energy density 
[28]. There have been some studies about the influence of the smoothness parameter a in the behavior 
of cU(z) [28], [30], [32]. In order to obtain the Dyer-Roeder like equation in f(R) gravity we follow [28]. 
First, we note that the terms containing the derivatives of in equation (4.54) comes from the 

transformation 4- — > 4- and the term with only c/^ comes from the Ricci focusing (4.24). Then 
following [29], [30] we introduce a mass-fraction a (smoothness parameter) of matter in the universe, and 
then we replace only in the Ricci focusing p — > ap. Thus following our arguments to obtain the equation 
(4.40), and considering the case Q r o =0we get 

d 2 4 (fl) - DR 4n m0 (i + zf + 3f(R)n k0 (i + zf + m DE *Qf> dd fw, dr 

(1 + z) dz* +(l + z > (a mQ (i + zr + f'(R)n k0 (i + z)2 + n DE ) dz 

, M(z)n m0 (i + z) 3 /(Jl)i dr. = „ , g 5] 

2(n m0 (i + z)3 + r(R)n kQ (i + z y + n DE ) A y '> 

where we denote the Dyer-Roeder distance in f(R) gravity by d^ R '' DR . This Dyer-Roeder also satis- 
fies the conditions (4.55) and (4.56), and it reduces to the standard form of GR in the case f(R) = R—2A. 

Is important to notice that the f(R) function and it's derivatives appear explicitly in (4.54) and (5.5), 
being this a big difference respect to the case of GR. Another feature for these theories is the role of the 
dark energy parameter £Ide, which contains the higher order terms of the field equations 4 . However, in 
order to show the effects of the higher order terms, is necessary, as we stressed above, fix either H(z) of 
f(R). Further studies are necessary to quantify these effects. 



6 Conclusions and Discussion 

We have studied the Geodesic Deviation Equation (GDE) in the metric formalism of f(R) gravity, start- 
ing for this from the form of the Ricci tensor and the Ricci scalar with the modified field equations. A 
generalized expression for the force term (Pirani equation) is obtained for f(R) gravity (for the FLRW 
universe), as well as the generalized GDE, both expressions reduce to the standard forms in the case 
f(R) = R — 2A. Then, we consider two particular cases, the GDE for fundamental observers and the 
GDE for null vector field past directed. 

The case of fundamental observers give us the Raychaudhuri equation and from this the standard form 
of the modified field equations with the FLRW metric as background. The null vector field case give 
us more important results: a generalization of the Mattig relation, the differential equation for the di- 
ametral angular distance in f(R) gravity, and as extension we consider also the Dyer-Roeder like equation. 

Another important result comes from the definition of the Ricci focusing: in the case of GR (in the 
FLRW space-time) we have a focusing for past-directed null geodesies for the condition n(p + p) > 0, 
here the condition is given by n(p + p)/f'(R) > 0. 

In the other hand, the general equation (4.3) is the starting point for important issues in f(R) grav- 
ity, as singularity theorems [33] and gravitational lensing [28]. We could study also the GDE for more 

4 From equation (4.44) the terms containing time derivatives of the Ricci scalar comes from the operator T>oq, equation 
(4.42), which gives the fourth order terms in the field equations of f(R) gravity. 
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general space-times (like Bianchi models) and for different forms of the energy-momentum tensor, in the 
context of f(R) gravity. 

Finally, the solution of equations (4.54) and (5.5) allows a different test for cosmological observables 
(areas, volumes, proper distances, etc.) in f(R) gravity, being an alternative method to constrain some 
models through cosmological probes. 
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